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O ' Abstract 



We derive a uniform approximation for semiclassical contributions of periodic orbits to 
the spectral density which is vahd for generic period-quadruphng bifurcations in systems with 
[ a mixed phase space. These bifurcations involve three periodic orbits which coalesce at the 
O . bifurcation. In the vicinity of the bifurcation the three orbits give a collective contribution to 
fH I the spectral density while the individual contributions of Gutzwiller's type would diverge at the 
O ■ bifurcation. The uniform approximation is obtained by mapping the action function onto the 
K*" ■ normal form corresponding to the bifurcation. The present article is a continuation of previous 
^ , work in which uniform approximations for generic period-m-tupling bifurcations with m 7^ 4 
H I were derived. 
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1 Introduction 



Semiclassical approximations for the density of states of a quantum system can be expressed in 
terms of classical periodic orbits |]T], 0, ^, ^, The form in which the periodic orbits contribute 
in these approximations is not unique; it rather depends on the characteristics of the classical 
motion. Up to now, complete approximations in terms of periodic orbits have been derived 
only for the cases of either integrable or globally chaotic classical motion. In the more general 
situation of a mixed phase space difficulties arise due to bifurcations, i. e. to the coalescence 
of two or more periodic orbits as the energy or some external parameter is varied. In the 
present article we describe the treatment of a generic type of bifurcation. We derive a uniform 
approximation for the joint contribution of orbits that participate in a period-quadrupling 
bifurcation. 

One way to derive semiclassical approximations is by starting from Feynman's path integral 
and evaluating all integrals semiclassically. The periodic-orbit contributions then arise from the 
stationary points of oscillatory integrals. If the periodic orbits are isolated the integrals can be 
evaluated by a stationary-phase approximation and one obtains Gutzwiller's trace formula. At 
bifurcations of periodic orbits, however, different stationary points coalesce and the stationary- 
phase approximation breaks down. This is because the different stationary points cannot be 
treated separately in the vicinity of the bifurcation. Instead one has to treat them collectively. 
This is done in terms of certain canonical integrals with the same structure of stationary points. 

The characteristic arrangement of stationary points in the vicinity of a bifurcation is de- 
scribed by its normal form. The generic normal forms for bifurcations in autonomous systems 
with two degrees of freedom or, equivalently, two-dimensional area-preserving maps have been 
classified by Meyer and Bruno |P, |^, H . They depend on the ratio m of the primitive periods of 
the periodic orbits which coalesce at the bifurcation. The corresponding bifurcations are named 
period-m-tupling bifurcations. Ozorio de Almeida and Hannay derived transitional (or local) 
approximations for the contributions of periodic orbits near generic bifurcations 0. These 
approximations are expressed in terms of canonical catastrophe diffraction integrals. They are 
transitional approximations because they are valid in the vicinity of the bifurcation. In farther 
distance to the bifurcation, however, they do not yield the correct amphtudes in Gutzwiller's 
approximation for the contribution of isolated periodic orbits. 

In two previous articles we extended the results of Ozorio de Almeida and Hannay by 
deriving uniform approximations which interpolate between the transitional approximation 
at the bifurcation and Gutzwiller's approximation for isolated periodic orbits [^, |Tl[ 



see 



also [Q). These uniform approximations were obtained by including higher-order corrections 
to the normal form expansion and then simplifying the integrals by appropriate coordinate 
transformations. The derivations were done for the cases m > 4 and m < 4, respectively. In 
the present paper we treat the remaining case m = 4. This case is more complicated than the 
other cases since it involves three periodic orbits whose action differences all increase with the 
same power of the parameter that describes the distance to the bifurcation. We use a different 
method for the derivation of the uniform approximation than previously: We apply techniques 
of catastrophe theory for obtaining uniform approximations for oscillatory integrals with almost 



coincident stationary points by performing a mapping onto the normal form |13, 14, 15]. The 



motivation is to provide a more elegant derivation; the results are, of course, the same for both 
approaches. 

In the following section we present the uniform approximation that is derived in appendix ^ 
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and discuss several limiting cases. The result is given for autonomous systems with two degrees 
of freedom and for two-dimensional area-preserving maps. We apply the uniform approximation 
numerically to the kicked top and discuss the limits of its validity. 



2 The uniform approximation 

The semiclassical contributions of periodic orbits to the spectral density can be obtained by 
expressing the density in terms of the trace of the (retarded) Green function, 



d{E)=Y^5{E-Er, 



-lmTrG'(E) , 



(1) 



and evaluating the trace semiclassically in the vicinity of the orbits. If the Green function is 
expressed in a mixed coordinate-momentum representation and the integrals over the compo- 
nents of the coordinates along the periodic orbits are carried out one arrives at an integral 
expression of the form 



d^{E) 



1 



Re 



dq' 



^ 1 dS 
J^rdE 



dpdq' 



exp 



-S{q,p,E)--qp--u^ 



(2) 



For a more detailed derivation of this integral see [T^, 11]. In (|^) the origin of the coordinate 
system is located on a central periodic orbit with repetition number r. Furthermore, p and q' 
are coordinates in a Poincare surface of section perpendicular to the orbit, and S{q',p,E) is 
the generating function for the r-th iterate of the Poincare map. It obeys the conditions 



OS 
dq' 



P 



dS 
dp 



dS 
'dE 



T 



(3) 



where the primed quantities are the final coordinates, the unprimed quantities the initial co- 
ordinates, and T is the time from initial to final point. The index ^ of the spectral density 
denotes the contributions from a group of orbits in the vicinity of the central periodic orbit. 
The periodic orbits are the solutions of 



dS 
dq' 



dS_ 

^ ' dp 



(4) 



and correspond to stationary points of the integral (|[). 

If the integral in (|^) is evaluated in stationary-phase approximation one obtains Gutzwiller's 
contributions of isolated periodic orbits. For a periodic orbit labelled by 7 this contribution is 
given by 



d^{E) 



A^{E) S^{E) 



cos 




where 



ME) 



T,{E) 



TrM^-2| 



(5) 



(6) 



Here S*^, T^, r^, and are, respectively, the action, period, repetition number, stability 
matrix, and Maslov index of the orbit. 
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In the vicinity of a bifurcation a stationary-phase evaluation of the integrals in (|^) is not 
appropriate. Instead one has to integrate collectively over all stationary points which are 
involved in the bifurcation. This is done by inserting the normal form of the generating function 
S{q',p, E) for the considered bifurcation into (||). For a generic period-quadrupling bifurcation 
the repetition number r is a multiple of 4; we denote in the following / = r/4. The normal form 
for this case is given by 

S{q\p, E) = So{E) + q'p - '-{f + - ^(g'^ + 2p^g'2 + - ^(g'^ - Gp^g'^ + /) 

= Sq{E) + g'p -el- ai^ - bl'^ cos(4<l>) , (7) 



where p = v2/cos$ and q' = v2/sin$. The parameter e is zero at the bifurcation. The 
normal form (|^ is obtained from an expansion of the Hamiltonian in the vicinity of the central 
periodic orbit, and 5*0 is the action of this orbit. 

The generic period-quadrupling bifurcation that is described by the normal form (|^) involves 
three periodic orbits, a central orbit and two satellite orbits. The bifurcation occurs in two 
different forms depending on the relative magnitude of the two coefficients a and 6 in (|^. In the 
case \a\ < \b\ there are two real orbits and one complex orbit before and after the bifurcation, the 
stable central orbit, an unstable satellite orbit and a complex satellite orbit. At the bifurcation 
one of the satellites becomes complex and the other one becomes real. For \a\ > \b\ the two 
satellite orbits are both complex on one side of the bifurcation and both real on the other side 
of the bifurcation, where one of them is stable and the other is unstable. The central orbit is 
real and stable on both sides of the bifurcation. 

In the vicinity of the bifurcation, i. e. for sufficiently small values of e, the contribution of the 
orbits can be described by the transitional approximation of Ozorio de Almeida and Hannay. 
It is obtained by approximating the pre-exponential factor in (|^) by its value at the origin and 
evaluating the integral with the normal form (|^ for the action. This yields the semiclassical 
contribution in terms of the diffraction catastrophe integral for the catastrophe Xg. In farther 
distance from the bifurcation the transitional approximation splits up into a sum of separate 
contributions of Gutzwiller's type. However, in this limit the semiclassical amplitudes come out 
with fixed relationships which are, in general, not in accordance with the periods and stabilities 
of the orbits. In more detail, the approximation is good as long as the following relations 
between the monodromy matrices and the periods of the orbits hold, 

4A^2iTrMo + A^2oTrMi = 8AS21 + 2A^2o 

4A^i2TrMo + A5ioTrM2 = 8A5i2 + 2A5io (8) 
A^2o Tr Ml + ASio Tr M2 = 2A^2o + 2A^io , 

where ASij = {Si — Sj)/2 and Tq = Ti = T2. The index denotes the central orbit and the 
indices 1 and 2 the two satellite orbits. Only two of the three equations in are independent. 
The relations ® follow from the normal form [cf. equations (|22|) and (0) in appendix^. 



For a general system they are only valid in the vicinity of the bifurcation. With increasing e 
they lose their validity, and the transitional approximation gradually becomes inaccurate. In 
order to obtain a formula which uniformly interpolates over the region from the bifurcation up 
to regimes where Gutzwiller's approximation is valid (without restrictions on the semiclassical 
amplitudes) one has to consider two modifications. First, higher-order corrections to the normal 
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form (|^) cannot be neglected any more. However, one can apply a mapping which brings the 
exponent in (H) back to the normal form. Secondly, one has to take into account the differences 
of the values of the exponential prefactor in (0) at the different stationary points. These steps 
are carried out in appendix 0. In the following we discuss the uniform approximation which is 
obtained there. 

The uniform approximation for the joint semiclassical contribution of orbits which are in- 
volved in a generic period-quadrupling bifurcation is given by 



1 [bP\ 
d^{E) ^ ---2 Re / dJ [To + a J + a2p] Jo Ht 
4/7rAi JO V / 

where J denotes the Bessel function of the first kind and 
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Equation (P) with definitions ([To|) is invariant under exchange of the indices 1 and 2. The index 
and (Tg, the sign of e, can be determined from the Maslov indices of the periodic orbits. The 
Maslov index of an unstable real satellite orbit is always u, that of a stable real satellite orbit is 
always u — ag and that of the central orbit is u + ag- The actions of the real orbits are ordered 
in the same way as their Maslov indices, i.e. Si > Sj if and only if Ui > uj. If both satellite 
orbits are complex the sign of e is given by = sign(S'i — Sq) = sign(S'2 — Sq). 



The coefficients (|T0|) depend only on quantities which enter also the individual contribution 
(|^). As an important consequence, the joint contribution is still invariant under canonical 
transformations. The transitional approximation is obtained if one keeps only the first of the 
three pre-exponential terms in the integrand. 

There are some special cases of values of a and b for which the integral can be evaluated 
analytically. They are discussed in appendix Numerically useful expressions that can be 
applied for arbitrary values of a and b are given in appendix |C[ In the following we discuss 
different limits of the integral @. If the action differences ASij are large in comparison to h 
then a replacement of the Bessel function by its leading asymptotic term and a stationary-phase 
evaluation of the integral yields a sum over Gutzwiller contributions (J^) for the real satellite 
orbits. Complex satellites do not contribute since they cannot be reached by a steepest-descent 
deformation of the integration manifold. (A detailed study of contributions of complex orbits 
near bifurcations can be found in This is sensible since the complex satellites, though 

having complex coordinates, still have real actions and semiclassical amplitudes [cf. equations 
(|22|) and (p^ j, and their isolated contributions would not be exponentially suppressed with 
h ^ 0. (In this respect, the role of the complex orbits resembles that of the complex satellite 



in a period-doubling bifurcation |TT[| .) 

The contribution of the central orbit is given by the leading semiclassical contribution from 
the boundary of the integral at / = 0. 
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In the opposite limit e = 0, i. e. at the bifurcation, all action differences vanish, 
leading-order semiclassical contribution of equation (01) is then given by 



The 



Tn 



4l^27r^h^\b\ 
if \a\ < \b\, and by 



=Re 



K 
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2\b\ 



cos 



5*0 



TT 



TT 

4' 



(11) 



(12) 



if \a\ > \b\. K{z) denotes the complete elliptic integral of the first kind. Tq and 5*0 are, 
respectively, the period and action of the orbits at the bifurcation, a and b are the coefficients 
in the normal form (|^) for e = 0, and Z = r/4 is the repetition number of the satellite orbits. The 
contributions (|Tl]) and (|12D are by an order h'^^"^ larger than the contribution of an isolated 
period orbit, i.e. the singularity index of the bifurcation is 1/2. Although it is not written 
explicitly, all quantities in ( pH]) and ( [T^ ) depend on the integer /. In detail, Tq^i = ITq^i=i, 
So,i = lSo^i=i, vi = 1^1=1 and a/ = /a/=i and bi = lbi=i. It follows that the amplitude of the 
contribution at the bifurcation decreases like l'^^"^ with increasing /. However, it cannot be 
expected that this approximation is good for arbitrarily large /. For longer periodic orbits 
bifurcations tend to occur more frequently. It is expected that for larger / there are other 
bifurcations which interfere with the considered bifurcation. 

The formulas of this section can also be applied, with minor modifications, to two-dimensio- 
nal area-preserving maps whose time-evolution is governed by the Floquet operator F. These 
maps correspond to systems with one degree of freedom whose Hamiltonian operator is periodic 
in time, H{t + T) = H{t), and F = U{T) is the unitary time-evolution operator for one period. 
The Floquet operator has unimodular eigenvalues e~^'^' with phases (pi that are called quasi- 
energies. The quasi-energies can be determined from a knowledge of the traces of powers of 
F. 

For maps the trace formula approximates Tr semiclassically instead of the level density. 
Furthermore, one has to pay attention to the following differences in comparison to autonomous 
systems with two degrees of freedom: (i) The orbits which contribute are those with a fixed 
period n, not those with a given energy E; (ii) the primitive periods have to be expressed in 
units of T and thus are integer valued; (iii) the action is not the reduced energy-dependent one, 
but depends on time (that is, on the number n); (iv) instead of taking twice the real part, the 
full complex contribution has to be taken; (v) the results further differ by a factor 2'n'h. 

It follows that the contribution C^"^ of a period-quadrupling bifurcation to Tr F^ is given 

by 

'bp' 



1 r°° ( bl^ \ ( i 

- I d/ |„ + aj + J„ [- j exp |-(S„ - il 



ITT 



(13) 



where the quantities a, b, e, ai and 0^2 are determined by (0) with Tq = n and Ai = 
n/ (r j y^l Tr A'/j — 2 1) . Here and Mj are the repetition number and monodromy matrix of 
the orbit, respectively. 
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3 Numerical results 



We now test the uniform approximation numerically on the example of a periodically kicked top 
r?, |TB|, |T^] and compare the results to those for the transitional approximation (valid close to 



the bifurcation) and the Gutzwiller approximation (which treats the orbits as being isolated). 
The kicked top is a dynamical system that involves the angular-momentum operators Jx, Jy, 
Jz which satisfy the usual commutation relations [Jfc, J;] = iekimJm, where h is set to unity. 
The evolution of the system conserves the total angular momentum + Jy + = j{j + 1). 
This introduces the quantum number j which fixes the dimension 2j + 1 of the Hilbert space, 
j + 1/2 further plays the role of the inverse of Planck's constant, and the semiclassical limit is 
reached by j — > oo. After normalization of the angular- momentum vector the phase space of 
the classical system is revealed as the unit sphere. 

The specific top that is considered here is described by the Floquet operator 



exp 



- iPzJz 



2j + l 



exp [—ipyJy] exp 



h 

^ 2j -\- \ Wx^x 



(14) 



This describes rotations by angles Pi and nonlinear rotations (torsions) of strength ki. For the 
study of bifurcations we hold the values of the Pi fixed [p^ = 0.3, Py = 1.0, Pz = 0.8) and vary 
k = kz = lOkx as a control parameter. The classical counterpart of the system is integrable for 
k = and displays well developed chaos at k = 5. 

At k = the top describes a linear rotation. In this situation the system has only two 
periodic orbits, both of period one and positioned at the intersection of the rotation axis with 
the spherical phase space. They are called the trivial periodic orbits. As k is increased, new 
orbits show up in bifurcations. The first two period-quadrupling bifurcations are encountered 
at A; = A;(i) = 1.0055 . . . and = fc^^) = 1.1954 .... Both have one of the trivial orbits in their 
center. The next period- quadrupling happens at k = k^^^ = 3.0336 .... The central orbit of 
this bifurcation is born at a smaller value of A; = fc*-^-* = 2.4497 ... in a tangent bifurcation 
together with an unstable partner. The form of the three period-quadrupling bifurcations is of 
type \a\ > \b\, i.e. on one side of the bifurcation {k < k^^''^'^^) both satellite orbits are complex 
and on the other side (fc > fc^^'^'^^) both satellites are real. 

Contributions of periodic orbits engaged in period-quadrupling bifurcations first show up 
in the trace Tr since the repetition number of the central orbit has to be a multiple of four. 
In the following we describe different semiclassical contributions to TrF^. In general, orbits of 
primitive period one, two, and four enter this trace. Complex orbits with the same primitive 
periods can also contribute, however they must be reachable by a steepest-descent contour 
deformation. All these orbits can participate in bifurcations. One has to deal with tangent 
bifurcations of orbits of primitive period one, two, and four, period-doubling bifurcations with 
central orbits of primitive period one and two, and period-quadrupling bifurcations with central 
orbits of primitive period one. For < k ^ 2.5, however, only few orbits are relevant, namely, 
those which participate in the period-quadrupling bifurcations at k^^^ and k^'^\ and in the 
tangent bifurcation at k^'^\ 

In figure |I| the trace Tr F*^ for j = 4 is plotted against k. The exact result is compared to 
three semiclassical approximations which, respectively, treat the two period-quadrupling bifur- 
cations by the uniform approximation, the transitional approximation, and the approximation 
that considers the orbits as isolated. The tangent bifurcation is described in all three cases by 
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Figure 1: The trace TrF^ of the 4th power of the Floquet operator (0) as a function of the 
control parameter k for j = 4. The plots show the exact quantum result, the uniform ap- 
proximation, the transitional approximation, and the approximation which considers the orbits 
as isolated. One observes that the isolated approximation diverges at the period-quadrupling 
bifurcations while the transitional and the uniform approximation behave regularly there. How- 
ever, the transitional approximation starts to break down far to the right of the bifurcations. 



the uniform approximation of . One observes that the isolated-orbit approximation diverges 
at the period-quadrupling bifurcations whereas the transitional and the uniform approximation 
behave regularly there. The transitional approximation, however, loses accuracy for the largest 
/c-values in the displayed range, which is most clearly seen in the real part of the trace. 

In figure |^ we investigate the behaviour of the semiclassical approximations as the semiclas- 
sical limit is approached, i. e. we fix the parameter k at k = 1.5 and increase the value of j. For 
low values of j the isolated approximation shows large deviations while the transitional and the 
uniform approximation are accurate. For larger values of j the sum of isolated contributions 
gains validity, since the effective Planck's constant l/(j -f 1/2) becomes small in comparison to 
the action differences of the orbits. The transitional approximation on the other hand becomes 
slightly more inaccurate, since the error in the semiclassical amplitudes shows up more strongly 
when the orbits can be considered as isolated. Figure ||(d) shows the deviation |Tr F^^ — Tr F^^\ 
of the semiclassical traces Tr F^^ from the exact trace Tr F^^. This function reveals a crossover 
between the transitional and the isolated approximation and displays the superiority of the 
uniform approximation over the whole range of j: The uniform approximation is up to an 
order of magnitude more accurate than the other two approximations. 

Figure ^(c) illustrates an oscillatory behaviour of Tr F^ which semiclassically originates 
from the interference of the contributions from the two period-quadrupling bifurcations. These 
contributions can be separated by considering the function 

"1 Jmax 

T(-)^S) = -. _ . E e-^-^TrF'^(j), (15) 

J — Jmin 

which has peaks at the positions of the actions of the periodic orbits. In its essence this function 
is a Fourier coefficient of TrF" with respect to j. It allows one to study the contributions of 
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Figure 2: The trace TrF*^ as a function of j for k = 1.5. (a,b) The real part of the exact quan- 
tum result, the uniform approximation, the transitional approximation, and the approximation 
which considers the orbits as isolated, (c) Illustration of the oscillations of the exact trace, (d) 
The error of the semiclassical traces, measured by the absolute value of their deviation from 
the exact one. 
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Figure 3: Peak heights \T^^\Sci{k))\'^ at the values of the actions of three periodic orbits as a 
function of the control parameter k. As k goes through k^^^ = 1.0055 . . . the three orbits coalesce 
in a generic period-quadrupling bifurcation. The plots show the exact quantum result, the 
uniform approximation, the transitional approximation and the approximation which considers 
the orbits as isolated. They are evaluated with jmin = 1 and jmax = 64. 
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periodic-orbit clusters individually in case that the action differences of orbits from different 
clusters are sufficiently large. A convenient testing tool then is the study of the peak height 
T^"'\S) at the value of the classical action 5*01 of a given orbit in a cluster as a function of k or 
for different values of jmin and jmax in ([T5|) . 

In figure ^ we show a quantum-mechanical evaluation of \T^^\Sci)\'^ for the three orbits 
involved in the period-quadrupling at k^^^ as the parameter k is steered across the bifurcation. 
We use jmin = 1 and jmax = 64. The exact quantum-mechanical curves are compared to results 
for the uniform, the transitional, and the isolated approximation. The uniform approximation 
is excellent over the whole range of k. It can hardly be distinguished from the quantum 
result. The approximation in terms of isolated orbits on the other hand fails completely. It 
diverges at the bifurcation and only gains validity again for values of k where the amplitudes 
of the contributions are already quite small. The transitional approximation is good at the 
bifurcation; however, as k is increased a clear deviation from the quantum result can be seen 
for the satellite orbits. This result shows again that there is a region where the uniform 
approximation is essential, since both other approximations, the transitional approximation 
and the approximation in terms of isolated orbits fail or are inaccurate. 

One can also investigate the peak heights for a given k and increasing jmin and jmax- In 
figure § the result is shown as a function of jmin with jmax = jmin + 63. Again, k is set to 
1.5. Deviations are visible for the transitional and the isolated approximation. Once more a 
crossover in the accuracies of these approximations is observed. The uniform approximation is 
by far superior over the whole range of jmin and cannot be distinguished from the exact result 
in the plots for the peak heights. 

Our numerical studies reveal that the approximation with isolated contributions is valid 
only far away from bifurcations or for sufficiently small values of Planck's constant while the 
transitional approximation does not gain validity in the semiclassical limit since it doesn't 
involve the correct semiclassical amplitudes. The uniform approximation however gives reliable 
predictions both close to a bifurcation as well as far away from it, and its accuracy increases in 
the semiclassical limit. 



4 Conclusions 

In the present paper we derived a uniform approximation for the joint contribution of periodic 
orbits that are involved in a generic period-quadrupling bifurcation. Together with the results 



of 1 10, 11] this completes the uniform treatment of generic bifurcations in autonomous systems 
with two degrees of freedom and two-dimensional area-preserving maps. Generic here means 
that these bifurcations are typically encountered in mixed systems without symmetries as one 
parameter of the system is varied. They are also called bifurcations of codimension one since 
only one parameter needs to be changed in order to bring the participating orbits into coales- 
cence. In systems with symmetries there can be additional kinds of bifurcations. Some of them 
can be described by small modifications of the formulas for generic bifurcations [^. 



In the following we discuss the limits of validity of the uniform approximations and possible 
extensions. Although in a generic situation one does not encounter any other form of coales- 
cence of periodic orbits than those discussed above, there are still cases in which the uniform 
approximations for generic bifurcations have to be modified. The reason for this is that pe- 
riodic orbits can undergo several subsequent bifurcations. One observes for instance that the 



11 



15 
12 
9 
6 
3 


150 
120 
90 
60 
30 


150 
120 
90 
60 
30 




exact 
isolated 
transitional 
uniform 




1 1 1 

unstable satellite 


1 






1 1 1 


1 



1 1 1 
stable satellite 


1 








'\/ 


1 1 1 


1 



isolated 
transitional 
uniform 




20 40 60 80 100 
j min 



20 40 60 80 100 

j min 



Figure 4: Peak heights \T^^\Sci)\'^ at the values of the actions of three periodic orbits as a 
function of jmin with jmax = imin + 63 and k — 1.5. The plots on the left show the exact quan- 
tum result, the uniform approximation, the transitional approximation, and the approximation 
which considers the orbits as isolated. The plots on the right show the deviation between the 
semiclassical and the exact peak heights. 
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iteration of the map corresponding to a normal form describes more periodic orbits of longer 
periods and allows for additional bifurcations. Another source of additional periodic orbits and 
bifurcations are higher-order terms in extended normal forms. If the bifurcations occur rapidly 
one after the other they cannot be considered separately, and instead all participating orbits 
have to be treated collectively. The next step beyond the isolated treatment of a bifurcation is 
the collective treatment of two subsequent bifurcations. Often the two consecutive bifurcations 
can also be considered as being part of a bifurcation of codimension two, since the two bifur- 
cations can be brought into coalescence by varying a second parameter. A collective treatment 
of two subsequent bifurcations then is necessary if one is sufficiently close to a bifurcation of 
codimension two in parameter space. The methods for obtaining uniform approximations for 
these cases are, in principle, the same as for bifurcations of codimension one; however, the 



normal forms are more complicated. Examples for these normal forms are given in [^, |21[] 



and uniform approximations for bifurcations of codimension two can be found in |22, E3l. It is 



even possible to completely semiclassically quantize a kicked top with low- dimensional Hilbert 



space by including uniform approximations for bifurcations of codimension two |24 



As longer and longer periodic orbits are considered, bifurcations tend to occur more rapidly 
and it is expected that then bifurcations of even higher codimensions become important. This 
makes the semiclassical treatment of mixed systems more and more complicated. Applications 
of semiclassical and uniform approximations in mixed systems are therefore most useful in 
cases where mainly the shortest periodic orbits of a system are needed, for example if one is 
interested in long-range fluctuations in a spectrum or if the contributions of long periodic orbits 
are suppressed. Applications of this kind are discussed e.g. in [p5|. 
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A Derivation of the uniform approximation 



We follow in this section a method from catastrophe theory for obtaining a uniform approxi- 
mation for an oscillatory integral with nearly coincident stationary points. By this method the 
exponent of an oscillating integrand is mapped onto a standard normal form with the same 
structure of stationary points. A description of the method and references to previous work 

A rigorous treatment including higher-order correction terms is given 



can be found in |Tj, |15 
in m. 



The contributions of the periodic orbits to the level density are contained in the integral (@) 



Re 



/oo 
dp g{q',p) exp 
-oo 



(16) 



where 



= S{q\p,E) - q'p 



(17) 
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and 



dS 

9iq,p) = 



BE \ dpdq' ■ ^^^^ 



The absolute sign of the mixed derivative of S inside the square root has been dropped since it 
is positive near the bifurcation (as can be seen from the normal form (|^)). It changes its sign 
only at conjugate points where the index u changes as well. By writing it without absolute sign 
the index u can be kept constant and the additional phase arises instead from the square root 
when its argument becomes negative. The energy dependence of the functions / and g is not 
written explicitly. 



We consider in the following the contributions of three periodic orbits to (|T6D that undergo a 
generic period-quadrupling bifurcation as the energy or an external parameter of the system is 
varied. We assume that any other stationary points of the exponent in ([161) which correspond 
to different periodic orbits are well separated from the stationary points that correspond to 
the periodic orbits which participate in the bifurcation. This means, for example, that the 
energy or other parameters of the system have to be limited to ranges in which the orbits do 
not participate in any further bifurcation. Under these conditions a uniform approximation is 
derived for the joint contribution of the three periodic orbits. 



Near the bifurcation the generating function S{q',p,E) in (|T^ is approximately given by 
the normal form (|^). From this normal form one can obtain properties of the periodic orbits 
as is done in the following. First the stationary points of (|1^) have to be determined. There is 
one stationary point at the origin which corresponds to the central periodic orbit. The other 
stationary points are determined conveniently in terms of canonical polar coordinates /, $ with 
p = -\/27cos$, q' = -\/27sin$. In terms of these coordinates the normal form is given by 

f{q'{I, ^),p{I, $)) = So-eI- al^ - bl^ cos(4$) , (19) 

and the stationary points of / are determined by the equations 

= sin(4$) , = -e-2[a + 6cos(4<l>)]/ . (20) 

There are altogether eight solutions for (ppj), four with cos(4$) = 1 corresponding to the satellite 
orbit which is labelled by 1 in the following, and four with cos(4$) = —1 corresponding to the 
orbit labelled by 2. The values of / at the stationary points follow as 

= ~2{a + a,,,b) ' ^^^^ 

where ai = 1 and a2 = —1. The satellite orbits are real if is positive, i. e. for o"e = —era where 
Ci = a + (Jib with i G {1, 2}, and we abbreviate the sign of a quantity x by ax- For negative 
the coordinates p and q become complex. The evaluation of (|T9|) at the stationary points leads 
to the values of the actions of the two satellite periodic orbits 

Si,2 = So+ f . (22) 
4(a + (Ti,2o) 

One can see from (|2^) that the action difference between any two of the periodic orbits increases 
proportional to for small e. The traces of the stability matrices can be determined from 

-1 , 



i dp dq' j \ ~^ dp dq' dp dq' dp^ dq'^ / ' ^ ^ 



d^S d^S d^Sd^S' 
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which has to be evaluated at the stationary points. It yields 



Tr Mo = 2 - Tr Mi = 2 + -^e^ , Tr M2 = 2 - -^^^ . (24) 

a + b a — b 



Both the actions and the traces of the stability matrices (|2^) are real quantities even for 
orbits with complex coordinates. 

Finally the Maslov indices of the orbits are determined. They are given by z/ + (ri„ — np)/2 
where ?t,„ and Up are the number of negative and positive eigenvalues of the matrix of second 
derivatives of /, respectively. They follow as 

uo = u + ae , 1^1,2 = ^(f^ci,2 - <yi,2<Jb) ■ (25) 



From (pi) it follows that for an unstable satellite orbit o"i^20'f) = "^ci 2 > whereas for a stable 



satellite orbit o"i^20"6 = ~<7ci,2- Using further the condition for real orbits crci,2 = ~'^£ o^i^ 
obtains 

Uq = U + ae , Us = J^-<ye, l^u = ^ , (26) 

where the indices u and s denote an unstable and stable real satellite orbit, respectively. 

We continue now with the evaluation of the uniform approximation. The equations (|2^) 
and (^) entail the conditions (^ which hold for e — > 0. If the distance to the bifurcation is 
increased (by changing the energy or a parameter of the system) higher-order corrections to 
the normal form cannot be neglected any more, and the dependence of S{q',p,E) on q' and 
p becomes more complicated. The main step in the derivation of the uniform approximation 
consists in the application of a coordinate transformation which then brings the exponent in 
([T6|) again into the form of the normal form (inside a region in which the stationary points are 
located) 

f{q',p) = F{Q',P) , (27) 

with 

F(Q', P) = So- ^{Q" + P') - ^{Q" + 2P'Q" + P') - ^-{Q" - 6P'Q" + P^) . (28) 

The parameters of F{Q', P) are chosen in such a way that the mapping from {q',p) to {Q', P) 
is one-to-one in a neighbourhood containing the stationary points. This can be achieved by 



mapping the stationary points [q'i.Pi) of the left-hand-side of ([27|) onto the stationary points 



{Q'l^, Pi) of the right-hand-side of (pTl) , which leads to the following condition 

f{q-,Pi) = F{Q',,P^) , (29) 

from which the parameters a and b can be determined. 

Condition (E^ is already fulfilled for the stationary point at the origin which corresponds 



to the central orbit. For the other stationary points the evaluation of (|29D leads back to (p^) 
where the parameters e, a and b now carry a tilde. Solving these equations for a and b results 
in 



16A5ioA52o ' 16A5ioA^2o 



(30) 
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where ASij = {Si — Sj)/2. Note that the third parameter e of the mapping is not fixed. 



The 



reason for this is that the form (^) contains one more parameter than is actually needed. By a 
simple scaling transformation Q' AQ', P AP one can change one of the three parameters 
into ±1. In analogy to (0) we define e by 



6^ 



2 - Tr M, 







(31) 



which measures the distance to the bifurcation. The sign of e is the same as the sign of e and 
can be determined from the Maslov indices (|26| ) of the central orbit and one real satellite orbit 
(if both satellite orbits are complex then = cr^). The difference between the quantities with 
tilde and without tilde is that those without tilde are obtained from a Taylor expansion around 
the central orbit whereas the quantities with tilde follow from the mapping. For e the 
mapping (^) approaches the identity transformation {Q',P) = {q',p) and the quantities with 
tilde approach the ones without tilde. 

The mapping (pTj) transforms the integral (p6|) into 



Re 



/OO /"OO 
dQ' / dP G{Q\ P) exp 
-OO J —OO 



where 



G{Q\P)=g{q\p) det 



d{q',p) 



(32) 



(33) 



and the determinant in (|33| ) is the Jacobian of the transformation which will be denoted by 
J{Q', P) in the following. 

The uniform approximation is obtained by writing the function G{Q', P) in the following 
form 

fjfP f)ZP f)TP f)p f)p 

GiQ', P) = ao - ai— - - + H,iQ', P)gQj + H,{Q\ P)— . (34) 

In order for this representation to be correct, the constants ao, . . . , as have to be determined 
such that the right-hand-side of (0) has the correct values at the stationary points. The last 
two terms in ( P^ vanish at the stationary points. They can be neglected since after inserting 
(|3^) into (|32|) they lead to terms which are of order h smaller than the other terms as can be 
seen by an integration by parts. Furthermore, the constant 0:3 can be set equal to zero, since 
after inserting ( p^ into (|32| ) the integral proportional to can be expressed in terms of the 
integrals that are proportional to the other as can be shown by another integration by parts. 

The remaining parameters ao,i,2 in (|3^) are obtained by an evaluation of (p3D at the station- 
ary points. For that purpose the value of the Jacobian J{Q' , P) at the stationary points has to 
be determined. This is done by differentiating ( [27| ) twice. With the notation {zi,Z2) := {q',p) 
and {Zi, Z2) := {Q', P) this results in 



dZ.dZ, 



E 



^ dZk dZi dZradZn 



1 -1i 



(35) 



from which the Jacobian J{Q', P) follows as 



J{Q',P) 



\ 



det 

det 



_dzkdzi_ 



(36) 
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By using (|23| ) this result can be expressed in the form 



J(Q',P) 



dpdq' 



Tr Mi - 2 



(37) 



where we defined 



Tr Mo = 2 - e 



;r2 



TrMi= 2 + ^52 
a + h 



TrMo 



86 . 



— £ 



(38) 



As before the quantities with tilde approach the ones without tilde as e — > 0. We continue now 
with the determination of the aj. By evaluating ( p3D at the stationary points with (|18|), (|3^ ) 
and (|37D one obtains the following conditions 



Tr M, - 2 
^ Tr Mi - 2 



(39) 



where /i and I2 are given by (pTD if the parameters on the right-hand-side of (^T]) are given a 
tilde, and Jq = 0. The constants Oq, ai and 02 follow from (^) and the definition (^T|) as 



ai 



0L1 



To 

4a (a + 6)2Ti 
— «o - - 



(40) 



8bP 



be \ (a + 6)(TrMi -2) 



+ 



be 



-8bP 



4(a 



;;2 



£2 



(S + 6) Ti 



86e 



=2 



+ 



^ (a-6)(TrM2-2 
(a - 6) T2 



2b \ (a + 6)(TrMi -2) 26 \ (5 - 6)(Tr M2 - 2 



86e 



~2 



This completely specifies the uniform approximation. The relations (|T0|) for the follow from 
(|40|) by using (^). The integral representation for the uniform approximation is obtained 
by changing the integration variables in (^) to canonical polar coordinates /, $ with P = 
\/27cos$, Q' = v^2Jsin$ and performing the integration over $. With the approximation for 
G{Q',P) that is discussed after (^) this results in the final expression (y). 



B Special cases 

B.l The case = (^i + S2)/2 

In case that the action of the central orbit is the mean of the actions of the satellite orbits, the 
diffraction integral in (^) can be evaluated analytically since a = 0. It is a special case of the 
form of the bifurcation with \a\ < \b\. The first satellite orbit is real if = — cr^ and the second 
satellite if = (Ti- We define in the following AS = (5*1 — 5'o)/2. The integral in (^ can be 
evaluated with the relation 

f d/ J.(/^) e.p(-,,/) = (e-K«)|) ff<;l (e"'4) , (41) 
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and the first two derivatives of tliis equation witli respect to rj. As before, cr^ denotes tlie sign 
of rj. Altogether one obtains the following result for the uniform approximation 



where 



\AS\ 



Re 



, 2 4v^ 4v^, 



Bi - icr^ 



An 



Bo 



+ 



'Ao _ _A^ _ A2 ' 
2 4^2 4^2, 



exp ( -So - y(z/ + ae 



Bi 

B2 
B3 



h\}1{z) H^yl{z) + {I ~ ia,)H\/, 
r(-) 



z)H\ 



/4 ( 



H 



1/4' 
(+) 



H[J>iz) H'S'iz) + H%iz) Hijliz) + 2(1 - ta,)H[j:^' (z) H^l 



-3/4 



'z)H 



-3/4^ 
(-) 



-3/4 



/4' 



1/4 



-3/4 I 



-3/4 V 



:i - i(T.)i7 



-3/4 > 



-3/4 ^ 



(42) 



(43) 



and 2 = \AS\/h. We have chosen the following notation for the Hankel functions: Hj^^^z) = 
Hl^\z) and H(-\z) = H(^\z). 

Near the bifurcation the actions of the two satellite orbits and the traces of the monodromy 
matrices are given by 



Si = So + 



4b 



S2 — So — 



4b 



and in the limit e — the contribution 



Tn 



16/7r2;i3/2|^|i/2 V4 
is obtained, which agrees with (|l^) for a = 0. 



TrMi,2 = 2 + 8e^ , 



Sn VT 



r^(.)cos(---. 



(44) 



(45) 



B.2 The case = (^o + S2)/2 



This is again a case where the integral in (^ can be evaluated analytically. This can be seen by 
considering the previous expression for the uniform approximation (|32D where the integration 
is performed in terms of the P and Q' variables. The constants a and b now satisfy the relation 
a = 3b, and the part of the double integral in (|32D with constant pre-exponential factor splits 
into a product of single integrals that can be evaluated. The bifurcation is an example for the 
case \a\ > \b\. Both satellite orbits are real if = — cr^, and complex if = cr^, and we define 
again AS = {Si — Sq)/2. The part of the integral in ( |5^ ) with constant pre-exponential factor 
is given by 



27r2| 


AS\Cf 









exp 



-2 AS] 



(46) 



where 



Ci — J-l/4 



AS 



h 



AS 



h 



^jo-ryr/S 



(47) 



and (Ti = sign (6) = sign (AS*). The terms with / or P in the exponential prefactor can be 



obtained from the first two derivatives of (|46D with respect to e. 



18 



We further define 



3/4 



h 



,-iCTg37r/8 



(48) 



and obtain for the uniform approximation (up to higher-order corrections in K) 



\AS\ 



Re 



'A 



Ai A 



Ao ^2 



+ 



An 



4 4v^ 



exp ( -S, - -u 



2 4 



(49) 



In the vicinity of the bifurcation the actions of the two satelhte orbits and the traces of the 
monodromy matrices are given by 



£2 £2 



Tr Ml = 2 + , Tr M2 = 2 - Ae^ , 



and in the hmit e — > the following contribution is obtained, 



32/7r2;i3/2|^|l/2 

This agrees with (O) for a = 3b. 



cos 



So IT 71 
— U 0"5 

h 2 4 



(50) 



(51) 



B.3 Further special cases 

We briefly discuss two additional cases that occur if 6 = or \b\ = \a\. 

In the case 6 = the actions of the two satellite orbits Si and S2 are identical. This can 
occur in integrable systems where both orbits are part of a torus. During the bifurcation this 
whole torus of orbits arises. With an integration by parts the formula (|^) reduces in the limit 
6 = to 

d^{E) ^ — -2Re / dJ [To + al] exp -(^o -el- ~al^) - -A , (52) 
4:LttTl Jo in 2 ) 

and thus can be expressed by a Fresnel integral. The constant a is given by 

a = lim (ai — ) • (53) 

6^0 V 2a J 

A discussion of semiclassical approximations for bifurcations in which a torus arises from a 
stable orbit is given in p6[| . 

The other case |6| = \d\ separates the two forms in which the generic period-quadrupling 
bifurcation can occur. The set of stationary points in the normal form (pH]) corresponding to 
one of the two satellite orbits goes to infinity as b approaches ±a. This normal form is not 
appropriate for a description of this case and correction terms have to be added to it. 
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C Series expansion of the uniform approximation 



The uniform approximation (P) consists of three integrals of type 



= / x" exp 

10 



—I'fx — lax 



Jo(a;^) da; 



(54) 



with V = 0, 1, 2. We now present power series in the coefficients a and 7 that are useful for a 
numerical evaluation of these integrals. In the case |a| > 1 (which corresponds to \h/a\ < 1 in 
the normal form) one expands the integrand around 7 = and uses the analytic continuation 
of the integral 6.621.1 in Expressing a hypergeometric function by its defining series and 
applying the duplication formula of the Gamma function one arrives at 



m=On=0 l^Q^J 



1 \ 2m r 
2^) 



■ - 



2m 



n\[m\\ 



(55) 



For \a\ < 1, i.e. \b/a\ > 1 in the normal form, one uses the integrals 6.699.1 and 6.699.2 in 



27| . Formally, a small imaginary part has to be added to the coefficients in order to assure 
convergence. At the end this imaginary part is sent to zero. After some transformations similar 
to those described above one arrives at 



11/ 00 



E 

n,m=0 



(V27A 



n^2m 



m\n\ 



(56) 



X 



sin (i±f±^7r) 



+ m 



r(i + m) 



la cos (i±f±^7r) 



2 / 



+ m 



r(| + m) 
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